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Summary 

A system of equations is developed for calculating 
spatial derivatives of flow quantities behind an expansion 
fan. For steady two-dimensional inviscid flow, equations 
in terms of the curvature of the streamline behind the fan 
are developed. Taylor series expansions of flow quan- 
tities within the fan are used and boundary conditions are 
satisfied to the first and second order so that the cur- 
vature of the characteristics in the fan may be determined. 
The system of linear equations for the spatial derivatives is 
then developed. 

An example of an application of the method is outlined 
for the problem of shock coalescence in which asym- 
metric effects are included as derivatives in the cir- 
cumferential direction. The solution of the coalescence 
problem may require values for spatial derivatives of the 
flow variables behind a resulting expansion fan. 

Introduction 

To solve supersonic- flow problems numerically, it 
often becomes necessary to know spatial derivatives of 
the flow quantities. Numerical computation of 
derivatives by use of difference formulas can sometimes 
cause difficulties or inaccuracies in the entire numerical 
scheme. The problem becomes less severe if analytical or 
semianalytical means of obtaining the derivatives can be 
found. Included in the first part of this paper is a method 
for finding spatial derivatives of flow quantities behind 
two-dimensional, nonuniform flow over a convex corner 
(i.e., an expansion corner). Taylor series expansions 
(ref. 1 ) are made for variables within the expansion fan. 
Boundary counditions are then used to obtain first- and 
second-order solutions for the variables and the cur- 
vature of the characteristics within the fan. Finally, equa- 
tions for the flow derivatives behind the fan are obtained 
in terms of the curvature of the streamline behind the fan. 

This method may be applied to the case of shock 
coalescence including asymmetric effects and can be 
combined with a sonic boom extrapolation program in 
reference 2. In that situation, coalescence of shocks often 
occurs at some distance from the body. This paper con- 
tains a description of the entire coalescence system in- 
cluding the resulting shocks, expansions, and contact sur- 
face. The asymmetric effects are felt through a second set 
of governing equations developed by taking the cross 
derivatives (derivatives in the cross plane) of the original 
governing equations. It is in solving the asymmetric part 
of the coalescence problem that it becomes necessary to 
know derivatives of the flow variables in the axial and 
radial directions. For calculations away from the body, 
some of the resulting shocks or expansions are often so 
weak that flow variables through them may be con- 
sidered continuous. In reference 2, a “continuous” ap- 
proximation for flow variables is made at the shock or 


the expansion of the opposite family which results from 
the solution of the coalescence problem. This approxima- 
tion limits the application of that method to cases for 
which the assumption of continuous flow variables is ap- 
propriate. The extension of that method to coalescence 
problems for which curved shocks of any strength may 
occur is discussed in reference 3, and the equations 
developed in the first part of this paper extend the ap- 
plicability of the method to coalescence problems in 
which expansion waves of finite strength occur. 

The second part of this paper describes the coales- 
cence problem and gives a description of how the equa- 
tions derived herein are incorporated into that solution. 

This work was done to satisfy in part the re- 
quirements for the degree of Doctor of Science at George 
Washington University, February 1983. 

Symbols 

A leading characteristic of expansion fan 
a speed of sound 

a( 0 ) zeroth-order term for speed of sound in region I 

B last characteristic of expansion fan 

C + any characteristic within expansion fan 

D vertex of convex (expansion) corner 

F\,F 2 , intersecting surfaces in shock coalescence 
F 3 .F 4 problem 

/] wall ahead of convex (expansion) corner 

f 2 wall behind convex (expansion) corner 

f' 2 slope of wall behind convex (expansion) corner 

f 2 change in slope in wall behind convex : 
(expansion) corner 

G = \p - g|(o) 

g | function defining leading characteristic of ex- 

pansion fan 

g 2 function defining trailing characteristic of ex- 

pansion fan 

g 2 curvature of trailing characteristic of expansion 

fan 

H enthalpy 

h slipstream surface in coalescence problem 

i unit vector in x direction 

j unit vector in y direction 

/ function defining streamline through expansion 

fan 



n direction normal to streamline 

0 order 

P dummy variables 

p pressure 

q m limit on maximum speed in inviscid, steady flow 
R gas constant 

r radial coordinate 

5 entropy 

T temperature 

u velocity in radial direction in polar coordinates 

u* velocity in x direction , in' Cartesian coordinates 

V velocity vector 

v velocity in angular direction in polar coordinates 

v* velocity in y direction in Cartesian coordinates 

Xi'i'), intersection of surfaces in coalescence problem 
R(¥) (see fig. 3) 

x,y Cartesian coordinates 

j3 shock angles; angle of trailing characteristic in 
fan 

7 ratio of specific heats 

6 flow angularity 

H Mach angle 

p density 

a unit vector in radial direction 

\j/ unit vector in angular direction 

a ,4/ polar coordinates 

V cross-flow direction in coalescence problem 
Subscripts: 

1 properties in region I; properties in region (1) in 
coalescence problem 

2 properties in region II; properties in region (2) in 
coalescence problem 

3 region ahead of expansion fan in coalescence 
problem 

4,5 regions behind shock and expansion fan in 
coalescence problem 

r first derivative with respect to r, d/dr 


rr second derivative with respect to r in coalescence 
problem, d^/dfi- 

x first derivative with respect to x, d/dx 

xx second derivative with respect to x in coalescence 

problem, d 2 /dx 2 

second derivative with respect to 'k in 
coalescence problem, dVd^ 

a first derivative with respect to a, d/da 

\p first derivative with respect to \p, d/d\p 

Superscripts: 

0 zeroth order 

1 first order, d/da 

A bar over a symbol indicates a quantity in region III 
(see fig. 1). A prime on a symbol indicates a first 
derivative. A double prime on a symbol indicates a sec- 
ond derivative. 

Part I — Development of Flow Derivatives 

Definition of Problem and Assumptions Made 

The expansion wave problem is illustrated in figure 1 . 
Assume the flow is around a convex corner with the 
vertex D at the origin (0,0). Let y = / t (x) for x<D (D~) 
and y=f 2 (x) for x > D (D + ) denote the sides of the 
corner. Flow ahead of the corner in region I is known to 
the characteristic line : C + at DA. The expansion fan 
(region III) is covered by C + lines from point D. The 
last C + line, DB, separates the fan from the flow field 
behind the corner. 



Figure 1. Flow over a convex (expansion) corner. 
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The problem is to determine the first derivatives of 
the flow quantities behind the expansion fan, that is, in 
region II at point D. In doing so, the first two terms in a 
series expansion in region II shall be defined along with 
the slope and curvature of the C + lines. The following 
conditions shall be used: 

1 . The flow quantities and their first derivatives (with 
respect to x and y) at point D (region I) ahead of the ex- 
pansion fan are used. These data should be consistent 
with the slope f[ (D~) and the curvature f\ '(D~) of the sur- 
face at point D ahead of the corner. 

2. The slope f ' 2 {D + ) and the curvature fi.(D + ) of 
the surface at point D behind the corner (region II) are 
used. 

3. All physical quantities such as velocity V, pressure 
p, and density p are continuous across the C + lines, but 
their normal derivatives may not be. 

4. In regions I and II, the solutions for the flow prop- 
erties are regular in x and y near the corner. In par- 
ticular, their partial derivatives d/dx and d/dy exist at 
point D. If P represents any flow quantity and j = 1 and 2 
for regions I and II, then the series expansion about (0,0) 
is 



Figure 2. Two-dimensional coordinate system. 


5. In region III, the flow quantities must change 
around point D, that is, they must be functions of ip as 
cr—0. With a bar over a symbol representing quantities in 
region III, 


Pj(x,y) = Pj(0,0 ) + x^Pj( 0,0) + y^Pji 0,0)- 


In polar coordinates ( o,ip ) (see fig. 2), this becomes 
Pj (o,iP) = Pj (0,0) + a ( cos ^ + sin * Pj (0,0) 


1 d -£ = o(- 

^ a dip \<7 


In region III, the power series in terms of o has the form 


Expansion (B) 


P(o,iP) = P<%) + ffPW(V') + 0.(<j2) (2) 


+ 0 ( ct 2 ) 

where O(a) indicates that the missing terms are of the 
order of a, that is, 

Expansion (A) 


where P0)( ip) = (d/do)P(0,ip). Expansions (A) and (B) 
are the two basic ansatzes needed for the analysis. They 
express the differences in the behavior of the solutions 
near D in regions I and II from those in region III, the ex- 
pansion fan. 


P; (ad) = P. (0) + aP. (1) ( 0,iA) + 0 (<j2) (1) 

J J J 

where Pj = Pj(0,0) and Pj(o,ip) is independent of as 
a— 0. 


Outline of solution. Basically, the procedure used in 
the analysis will be to substitute expansions (A) and (B) 
(eqs. (1) and (2)) into the governing equations and then to 
compare the coefficients of like powers of a. An outline 
of the solution procedure is the following: 


Stated alternatively, 


dx 


, d . . ,ia 

COS - sin V'— T-r 
oo o oy/ 


_3_ 

dy 


. , d , , 1 d 

sin \p — + cos - -r-r 
do o dip 


If d/dx and d/dy are finite as a— 0, it is necessary that 
(1/a) (d/d\ p) is finite as a— 0, that is, d/d\p~0 as a— 0. 


Step (1): Set up the general equations of two- 
dimensional motion with nonuniform enthalpy and 
entropy. 

Step (2): Determine the slope and curvature of the 
first characteristic line in region III (DA) and the bound- 
ary conditions between regions I and III along DA. This 
yields the first term in expansion (B). 

Step (3): Determine the flow properties at point D in 
the expansion fan, those behind the fan, and the slope of 
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the last C + line (DB) from the slope fi(.D+) of the sur- 
face^. 

Step (4): Determine the second term in expansion (B) 
(eq. (2)) in the expansion fan and the curvature of the 
characteristic lines. 

Step (5): Apply boundary conditions along DB (con- 
nection between regions II and III) and determine the 
first derivatives (with respect to x and _y) of the flow prop- 
erties in region II, the derivative being compatible with 
the curvature of the surface behind the fan. 

Basic definitions and relationships used in solution. 
Included below are a few definitions and relationships 
which are used in the subsequent derivations. The com- 
ponents of the velocity along the radial unit vector a and 
the circumferential unit vector \p are u and v. Thus, 

V = m( \p) + v\p(\p) 

The components u and v are in turn expressed as power 
series in a as follows: 

u = V* o(f) = u ( %) 4- au ( '\f) + 0(o2) 
and 

v = V . 4s(f) = v (0 W) + av (l) 0A) + 0 (<t2) 

In regions I and II, V(cr = 0) exists; hence, 
u { %) = V(<r = 0) • o(f) 
v (0) W = V(cr = 0) * <//('/') 
and 

u iV] (f) = V^O 
v {V> (f) = VM) • d(f) 

Also note that 

uj® = V(a = 0) • yf/(ip) = v (0) 
and 

vf = -V • o#) = 

Solution Procedure 

Step 1. In polar coordinates (o,f), let u be the radial 
velocity and v be the circumferential velocity. The govern- 
ing equations in regions I, II, and III are as follows: 

Continuity 

U a + IT + o" ~ p (j*Po + ~oP<p) = 0 (3) 


where p is the density. 
o-momentum 

V V2 1 

uu a "I" a u 'l / ~ ~a ~ ~ p (4a) 

where p is pressure. 
i p-momentum 

uv a + T v^+ — = -—p^ (4b) 

Energy 

uH a + 1^ = 0 (5) 

where H is the total enthalpy defined by 

„ a 2 , « 2 + v 2 

tl — 7 H ~ — - 

y - 1 2 

and a is the speed of sound. An additional equation, 
valid along streamlines, is 

uS a + -^S^ = 0 (6) 

where S is entropy. Eliminating p and p from the 


momentum and continuity equations yields 



1 / d d \/ u 2 + 

V 2 ^ 

(7) 

cm a + u + - a2 (™ da + v di/ )( 2 

and 



(u2 + v 2 ) (u^ - ov a - v) = uHj, - ovH a - T(uSj, - 

voS a ) 

(8) 


where T is the temperature. Equations (5), (6), (7), and 
(8) are the forms of the governing equations used in this 
derivation. 

Step 2. Let f = gi(cr), and thus G = \p - g\(o), be the 
equation for the slope of the first characteristic line (DA) 
in region III. The velocity perpendicular to DA is a, that 
is, 

y.ZC = bmM + (.,/ gl , ±c (9) 
I TC I [(l/olHgiWj 4 

where =dg\/dd. Since Vi < 0, the negative sign on a is 
chosen. Expanding both sides of equation (9), one gets 

[-^ 0) ag{(a) + vj + <7V ] (1) + . . .Jjl -jo2[g{(o)]2 
+ ...(= + C»(a 2 )] 
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But 


since 




dvi 

do 


= &?* = V0*i = 


d\l/ da 


'1^' 


-ufhi 




Therefore, 

(a) + Vl - a«f Vi]{l - ja 2 [gi(<j)] 2 | 

= -fa® + aa® + <3(ff2)J (10) 

Equating a for the zeroth order,, one gets 


implying that gj(0) = tan (6\ + m) at point D (i.e., it is a 
characteristic line). Equating a for the first order, one 
gets 

-2M®gi(0) = (0 = ^(0)) 


or 


g;(0) 


a ( J 1) g]( 0) 


2 u 


( 0 ) 


(ID 


Define 0j = gi(0) to be the inclination of DA at D. Then, 
the curvature of DA at D is gj(0) = 30/3ff| g Also 
note that 



da I (a,^i)' 

da 



da\{x,y) 

dx 


COS 01 


+ 


fai(x,y) 

dy 


sin 


at x = 0 and y = 0 in region 1. 

Along the boundary DA (0 = g\(a) = 0i + ffg|(0) + 
. . .) the flow quantities are continuous, that is, 


P 1 (di&lO 7 )) = P(o,g\(a)) 


or, from ansatz (A) and ansatz (B) (eqs. (1) and (2)), 
P\ + ^/W)) + 0(a2) = 

P (0 \g\(°)) + oP ( \«o)) + 0(o2) (12) 


at (7 = 0. Recall that in region I, all flow properties 
are finite and single valued as a — 0; however, in region 
III, the properties are not single valued but are functions 
of 0 at a = 0. For the velocity V, 

Vl((7,gi(a)) = V(a,gi(a)) 

The velocity at point D, as indicated in equation (13), is 

v; o) = v (o) (g,(o)) = v (o Vi) 

= «(°>M) a (0,) + v(°\h) 0(0i) 
and the velocity to the first order is 

V U) (0i) =[«JV(0) + h (1) (0,) - v (0) g;(0)] a(0,) 

+ [w (0) g'(0) + v^V(O) + v <l> (i/'i)]i/'(i/' 1 ) 
where the properties da/d 0 = 0,30/30 =-<r, da/da = 
0, and d\f//da = 0 are used. 


Step 3. From the application of expansion (B) 
(eq. (2)) to governing equations (3) to (6) and the 
equating of the coefficients of <r®, the resulting zeroth- 
order equations are 


o 

II 

(15) 

s< 0) =o 

(16) 


_(0) _(0) 2 _(0) _(0) _(0) (-(Of _(0) 2 U0) „ 

w ’ cr ’ -ir’v y ’u+\a - v J F ' = 0 (17) 

0 0 

(#>’ + pMViS/2 _ # i) , 5,0) H<0)_ 5,0) ,*0) M 

'\ .30 / f f 


Equations (15) and (16) imply that H and S are in- 
dependent of 0; thus, 


where P represents the general flow quantity. Therefore, 
Pi = P i0 \i i) (13) 


ii 

Sc 

(19) 

s, = s 2 

(20) 


and 


From equations (15) and (16), equation (18) becomes 


^ 1 1) (0 1 ) = ^° } (0l) *1(0) + (14) 


-(0) _(0) 
u - v 


= 0 


( 21 ) 
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From this, equation (17) may be written as 


From equations (24) and (25), 


(sW-p® 1 )^ # 0> ) = 0 


A - 1 A 


The zeroth-order term of equation (4b) gives 
v< 0| (v«+ B®>)= ‘ 

' $ ' p(°) W 

Since ^ 0 at point D in region III, then 

v< 0) + « (0) * 0 
v 

and equation (22) becomes 

5 W -V <0 >‘.0 


From the direction of the flow, the negative sign is 
chosen; thus, 

« (0) = (23) 

which, as expected, defines each local (a —0) radial line 
to be a characteristic. Equations (19) and (23) may be used 
to relate T7(°) to 3(°). From the definition of qfa as the 
limit on maximum speed in inviscid, steady flow, one 
obtains 

jjW + yw | 3 w jjW + ( 

2 7-1 2 y- 1 




u (0) n (0) /7 + IV/2 

Along a streamline, defined by a = one ob- 

tains (at the boundary DA) 

cti = l{^ = g\(ax),a x ) (29) 

Expanding from the boundary DA along the streamline 

aj = /0/sffi) = / (0) ty) + < r/W ($) + ... (30) 

one gets to the zeroth-order and first-order equations 

0 = / (0) Wq) 


1 = / (1) (^,) 

Therefore, for a streamline, 


IT = (Tj / (x) (^) + 0((j2) 


Recall that 


« = = £L ;(0 = J_ ,(i) 

v <7 dxp a yj/ 1 ( 1 ) \p 




Equation (21) then becomes 
H (0) = v< 0) = -a (0) = - 


and after integration 




u ( %) = q m smd> (25) 

where 4> = t( 7 -l)/( 7 + l)] l/j (a- \p) and 7 is the ratio of 
specific heats. From boundary condition (6), 

n (0) «7) = «®= q m sin (a - W 

and 

» ‘ ( sin ‘‘ + (26) 


1 c// {1) _ U (0) /y+ 1\'/ 2 , . 

/(» d* y (0) "(7-I/ ^ 


In/ 1 ') - in 


7 - 1 COS <^>1 


From the definition of cr given above, 


= /(1) ( ,/a = / cos <A - 
a, \cos <t> x J 


■(7+l)/(7-l) 


To determine the slope of the last characteristic in 
region III ( DB ), write the equation 

'P = = &(0) + °g' 2 (0) + 0(o2) 

and define £2(0) = ^2 as the inclination of the 



characteristic at D. The boundary conditions along DB 
between regions III and II are similar to those between 
regions I and III along DA, that is, the flow quantities 
along DB are continuous. These conditions imply the 
following: 


Similarly, from equation (6) and the boundary conditions 
at DA for S, 


S (1) (^,) = sj 0 


lim 

da 




D 2 (a,g2(a + )) = P{o,g 2 (o-)) 
Pl(a — 0) = P 2 = P (0) (h) 
u 2 (a = 0) = u 2 = 
v 2 ( ff = 0) = v 2 = v (0) ty- 2 ) 


cos + — l - sin \p, 
dx 1 


dS i 
dy 


x=0 

o 


(37) 


and 


s (1) W 

cd) 


x(T + ')/(7-l) , 

cos <j> \ 1 

COS 01/ 


l (!) w 


(38) 


An implicit equation for \ p 2 is derived using the slope of 
the wall behind point D, y = f 2 (x). The slope is 

f . . _ dy_ _ -u sin \p + v cos \p 
dx u cos ^ + v sin \p 

/'(0) = ,-»2li,n^2 + V2COS^ 2 (33) 

u 2 cos \p 2 + v 2 sin \p 2 

The results obtained in this step are identical to those 
in a centered expansion wave separating two uniform 
regions (ref. 4). They have been included so that the next 
step provides a continuation to take care of the non- 
uniformity of the flow fields ahead of and behind the 
expansion fan. 


The first-order equations for w (1) and i> (1) are derived 
from the coefficients of a in equations (7) and (8) The 
first-order equation from equation (7) is 

*<»(-#> « (0 >) + fP(2a^ - /» 2 - v (0) u<®) 

+ vj>(o< 0 > 2 - v<°> 2 ) + v<»(-u< 0 > v (0 > - u< 0) uf 
-2v (0 > + (39) 

The following relations, derived in step 3, are used to 
simplify equation (39) and a similar equation for equa- 
tion (8). From equation (21), 


Step 4. In this step, values for the next-order terms 
(w^O/,), v (1, 0/<), // (1) W0, and are determined. 

Also determined is the curvature of the characteristics at 
point D by comparing coefficients of a in equations (7), 
(8), (5), and (6). 



v<°) 


and from equation (22), 




From equation (5), one obtains 

u (0) (i/0 H ( %) + v< 0) H%) = 0 (34) 

4 

and therefore 



The boundary condition at DA yields 




lim dH\ 
ff -*° da 


4= 4, 


dH x , ^ d//] . , 

- cos h 1 sin \p l 


dx 


dy 


x=0 

y = o 


Because of flow direction, choose (eq. (23)) 

= -#> 

Likewise, equations (23) and (27) give 



-v<°) 


As given in equation (25), 


W’ = q m sin </> 

where <j> = [(7-1 )/(t+ 1)] 1/2 (a-\p) and d<p =-[( 7-1)/ 
(7+ l)]'/z d\p. Differentiating cft°\ one obtains 



1^1 w (0) 
7+ 1 


7 



and differentiating equation (23) gives 


v(°) 




Substituting for v^\ one then arrives at the following 
expression for in equation (39): 

+ v^ = _ 2 _S (0) 

4> 7+1 

From the definition of enthalpy and dH/da, one obtains 
a ([)2 = ( 7-1 )(/70 - u<°) - #) vO) (40) 

Solving equation (39) for v (l) yields 


v<» = - 


+ 


«J ) + t7 (1) {fl ( ° )2 -i7 (0)2 [(37-l)/(7+l)]} 


2 

7-1 H (X) 
7 + 1 y(0) 


20<°> v<°> 


(41) 


A similar procedure for equation ( 8 ) yields 


p(i)= 4 ' )_ i) (//(i) ' r(0)s(i)) (42) 


By equating equations (41) and (42), the following equa- 
tion for « (1) is obtained: 


v 


37-1 


2(7+1) \7 - 1 


=( 0 )^CX) 


2 \7 + 1 


' - 3 + r"’t 


f+l 2 #) 2 v< 0) 


The integrating factor for this equation is 


(43) 


9= exp 


37 ~ 
2(7 + 

V /2 


+ 


i_/t + i\ /2 

1) \7— 1/ • 


f cot <j> dp 


f tan tp dip 


T/V-^y 

2 \7+ 1/ 

—■ ) r S tan 4> d<p - 4- S cot tp dtp 
2(7-1) 2 J 


9= (cos <p) <3 > 1) / i2(T " 1)1 (sin0)’ 


Zi 


(44) 


Thus, 


m ( 1) = (cos (pj iy <7 '' (sin <P) Vl ( J' x dip 

\^i 

) 9 


+ / 


(*i) 


( 45 ) 


where 


5<%i) = u^iiPO-ufHh)?'^) 

+ v <0) ('/'i) g{(°) = w ' 1) (’Ai) (46) 


and 


(COS ^-'^-Dksin 0)'* ^ dip 
7+1 2v (u) 


Hi / 1 


,(7 + 


l)/(7-l)^ (COS 0)C 3 t+ 5 )/[2(>-O1 


Qm [ cos ^ (sin <f>) Vz 


d<t> 


„ * ,,-(37-11/12(7-1)1, . T (0) S (1) ,, 

f diP = f (cos <p) (sin 0) /; "„_ im - diP 


2v<0) 


1 „(cos<p) y + 


2y R 


If 7 = 7/5, then 


( 1+4 

\cos <Px) 


1/12(7-01 


f 


(sin <P) 


>/2 


- dtp 


/ \ / \6 

^ dp = 2 (Sj^co^r) ( sin ' /2 tp - sints 0i) (47) 


and 


^ .. 5S,<? W / 1 V / 2 


f tJ& 2 dp = 


14 i? \cos 4>i/ \ 7 


— cos2 tp sin-'+ <P 


8 2 

-I sin'+ 0 - — cos2 0i sin- +2 0 ( 

35 7 


sin '+ 0. 

35 1 


(48) 


Knowing n* 1 ), one can evaluate equation (43) for u ( J), 
equation (42) for and equation (40) for a*') 2 . 

To evaluate the curvature of the characteristic C + in 
region III at point D with inclination ip 0 , define the 
curvature ,with the equation ip = g(o,ipo), that is, p 0 = 
g(0,ipo). For the first and last C + lines in region III, 
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gl(a) = g(a,\l/[) 

gi(<y) = g(a,M 

Given this definition for g(a, xp 0 ), the dependence of g on 
i/'q shall not be shown. The normal component of veloc- 
ity is a(a,g(a)), or 

(ua + v^) » VhA + g(g)] = - 
|V[^-g(<r)]| 

Ug'ja) - (War) = _ 

[g'(a)2 + (1/ u 2 )] Zi 


or, equating a and xp components 

I/O) = v< 0) )g' + i/0) 

vG) = g' + y0) 

yields 

4 1) (^2) - S^to) (52) 

v<% 2 ) = ^ W 2 ) + (*f + t/ (0) )|^ 2 g'^2) (53) 

Note that 


Recall that xp = g(o) = g(0) + <rg(l)(0) + . . . , where g(0) 
= V'o is the angle at a = 0, or 


fp'Hh) 


lim 
o~~ 0 


dH 2 

da 


'l'='P 2 


dH 2 

dx 


COS xp 2 + 


dH 2 

dy 


sin p 2 


[Mg'(ff) - (v/o)]o[\ +<j 2 g 2( (J )]!/2 


Since we know the slope 0 2 of the streamline, we obtain 


= #># 0 ) + a[g\0)df + aO^o)] 

Expanding the left side and equating the coefficients of a, 
one obtains 


^°\h)g'( o) = + v^Wo) 

a (1) (i/> 0 ) + v (1) (^q) 
M (0) ( V'o) 


Curvature = g'(0) = 


Also note that 


«(D 2 = lim „ = ''m 2 o< 0 > ^ = 2a (0) «0>M 

d# a— u da 


Therefore, 




q0%) 


Step 5. The boundary conditions along Z)5, the last 
C + line in region III, are developed in the same manner 
as those along DA. Along DB, xp = g 2 (a) and g2(0) = 
xp 2 . The term g' 2 ( 0) is given by equation (49) with 1 p 0 
replaced by xp 2 . Applying the boundary conditions yields 

HfWi) = «2°W (50) 

$\h) = $iHh) ( 51 ) 

( I/O) a + yO) ^ = (u^ 0 + «*°) xp + pG) xp - v^cr^g' 
+ i/O) a + y(') \p 


Similarly, 


dH 2 . . dH 

- 1 -sin (09 - w 2 ) = — 

dn 

da 

dH 2 

/4 1) (^ 2 ) 

dn 

£) sin (0 2 + M 2 ) 

dS 2 1 

s ( 2 V 2 ) 


(54) 


dn 


D sin (0 2 + M 2 )- 


(55) 


Denote the velocity components in the x and y directions 
in Cartesian coordinates by u* and v*; then, 


«4 0 (* 2 ) = V 2 , <J i(a,^ 2 ) • °# 2 ) 


x = 0 

y = 0 


= cos xP 2 + u *, sin 
+ cos xp 2 + V* y sin xp-^j 


cos Xp 2 


x=0 

y = 0 


sin 1 p 2 


or 


u ( 2 ]) (xp 2 ) = u£, (0,0) cos 2 xp 2 

cos xp 2 sin xp 2 [>(0,0) + V^(0,0)] 


+ 


+ v* y (0,0) sin 2 xp 2 


(56) 
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Similarly, 

= sin \p 2 cos <A 2 [-<(0,0) + v<(0,0)] 

- <(0,0) sin2 \j /2 + < (0,0) cos2 i 2 (57) 


where 

V 2 (x,y) = V 2 (0,0) + a {[*<0,0) + ,<0,0)] i 
+ [<0,0) * + v*(0,0),]j| + 00/2) 


Equations (7) and (8) or (3) and (4) in Cartesian coor- 
dinates are equivalent to characteristic equations. Since 
the boundary conditions along DB are used, the equa- 
tions along the C + lines are fulfilled to the a order. Only 
one independent equation is available from equations (7) 
and (8) or (3) and (4). Choosing equation (4) we get 


-g 2 (v x -u y ) = 




| where q(d/dn) = u(d/dy ) - v(3/dx) and d/dn is the 
derivative normal to the streamline. The leading term 
of this equation (<j°) yields 


«72(0,0)[ M ;(0,0) -<0,0)] 

_ dH 2 (0,0) _ <0,0) dS 2 ( 0,0) 


dn 


y R dn 


(58) 


where R is the gas constant. 

One more equation, ,=/ 2 (x), is available from the 
boundary condition along the wall. Thus, 


m 2 / 2 'W - v 2 = 0 


d_ 

dx 


( 0 , 0 ) 


[«*/'(*) - V*] 


along , = f z (*) yields 

[<( 0 , 0 ) + <( 0 , 0 )/ 2 '( 0 )]/'( 0 ) 

' [<(°>0) + v <(0,0) /' (0)] = -k*(0,0)<0) (59) 


Equations (56), (57), (58), and (59) are four equations for 
<(°,0), <(0,0), <( 0,0), and <(0,0) in region II. In 
coefficient matrix form, these become 


COs2 \J/ 2 

sin \p 2 cos \p 2 

sin i p 2 cos \p 2 

sin 2 

< 

-sin \ j/ 2 cos \j/ 2 

-sin2^ 2 

cos2 \j/ 2 

sin xp 2 cos 1^2 

U 2y 

0 

1 

-1 

0 


4 

iff 

-1 


V 2y 


< ih) 

< (* 2 ) 

1 dH 2 A dSl 

Q 2 dn 7 R dn 

-k* 2 }" 2 ( 0) 

The determinate of the coefficient matrix yields 



cos2 1 p 

sin 2\p 

sin2^ 

det = 

-sin \p cos \j/ 

cos 2\l/ 

sin \p cos \l 


f 

J 2 

iff-\ 

-f' 

J 2 

-sin 2h 2 j 0 
2cos2 0 2 

(0 <u 2 < 90°) 


Part II — Application of the Method to the 
Case of Shock Coalescence Including 
Asymmetric Effects 

The coalescence system is briefly described herein. 
For a detailed derivation of the governing equations 
of the system, see reference 2. 

The first two shock surfaces, T 7 , (r,'k) and F 2 {r,^), 
coalesce to form a resultant shock surface F 3 (/\'k), a con- 
tact surface h{r ,¥), and a weak (isentropic) shock or ex- 
pansion of the opposite family F 4 (rji). The solution 
developed in part I is applicable to the case for which F 4 is 
an expansion wave, as shown in figure 3. 

The problem is first solved axisymmetrically, and at 
some distance away from the body, it may be treated as a 
two-dimensional problem, as shown in figure 4. All con- 
ditions in regions (1), (2), and (3) are assumed to be 
known, as is the point of coalescence. The unknowns for 
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(r) -Resultant shock 



h <x) - Slipstream 


Figure 4. Two-dimensional representation of shock coalescence. 



this problem are the following: velocities u 4 , v 4 , «*, and 
v*; pressures p 4 and p 5 ; temperatures T 4 and T 5 ; the 
angle of the resultant shock with respect to the upstream 
flow 183; and the angle of the last characteristic of the ex- 
pansion fan (84. The system of equations consists of the 
four shock equations across F 3 , the four isentropic condi- 
tions across the expansion wave F 4 , and the matching 
pressure and flow direction at the slipstream surface h. 

The system is closed, with 10 unknowns and 10 equa- 
tions. To derive the asymmetric system of governing 
equations, the system is treated as five intersecting sur- 
faces along A r ('^), R(ft). (See fig. 3.) The governing set of 
asymmetric equations are the 5 intersection equations, 
the 5 shock equations at F 3 (l continuity, 1 energy, and 3 
momentum), 5 equations across the expansion wave F 4 , 
and 3 slipstream equations at h for a total of 18 equa- 
tions. The second derivative with respect to 'k (the 
circumferential direction) is taken for each of the govern- 
ing equations along the corresponding surface. The equa- 
tions are then reduced to the 'k = 0 plane, resulting in 18 
equations and 18 unknowns of the form m*¥|* = 0 > 
,P^| +=0 . an d so forth. The asymmetric system would be 
closed except that in the asymmetric equations there ap- 
pear spatial derivatives of u*, v*, p, and T in both 
regions (4) and (5). The method derived in part I is now 
applied to provide a method for obtaining the spatial 
derivatives. Equations (56), (57), (58), and (59) give u 4x , 
v| x > u 4r , and v 4r in terms of the curvature of the 
streamline behind the expansion fan ( h ** for this problem). 

The Euler equations in the plane of symmetry (ty = 0) 
are valid in regions (4) and (5). The following four equa- 
tions are applied in region (5). Because the curvature of 
the characteristics in the expansion fan have already been 
determined by the method developed in part I, only equa- 
tions (60) to (62) need to be applied in region (4). 

x -momentum 

p$u*u* x + p$v*u* r + RT 5 p 5x = 0 (60) 

r-momentum 

P5 u * v *x + P5 v 5 v *r = - 8 P 5 ~ (61) 

Continuity 

P 5 T 5 u* 5x + Ps T 5 ,v* 5r - u$p 5 T 5x - v 5 */? 5 T 5r 

+ + v*sT s p 5r = - - P - ~ v* (62) 

Energy 

(Us) 2 T 5 U*s x + U*V*V* X + U*V*U* r + (v*) 2 Vv 

+ ~l U * T 5X + ~f* T Sr = gy* (63) 


The tangential derivative of the three properties along 
the slipstream h yield 

PAx + PAr^x = P5x + Psr^x (64) 


^Ax^x + ua^xx ~ VAx + u 4r h 2 -v 4r h x — 0 (65) 

u 5x^x "k w 5 ^xx ~ v Sx + u 5r^x ~ v 5 r^x = 0 (66) 


The remaining equations at / 3 are derived from the 
following four shock equations: 

Continuity of mass 

P\TsQ\ = PsT\Q 5 


r-momentum 


(«5 - «i) Fir + (v? - vft = 0 


Energy 

(U* + Mj) (1/5 - M*) + (v* + V*) (V* - V] 1 ) 

, 2yR(T 5 -T{) _ Q 
y- 1 


Combination of all shock equations 


■(’ ‘ #) - 


2yTA\+F] + -f = (7+DQ1Q3 - (7-1)0? 


where Q = u*-v*F r Taking the tangential derivative of 
each of these equations along the resultant shock F 3 
yields the following four equations: 


(P\Ts x Q\ - Ps x T\Qs - PsT\U* x + P5TiV*xF 3r )F ir 

+ P\T$ r Ql ~ P5rT\Qs - p$T\U\ r 

+ Ps T\ V 3r F 3r + PsT\Fs rr + P\ TsV*F 3rr 

= (P5 T lxQ5 ~ P\x T sQ\ ~ P\TiQ\x)Fsr 
+ Ps T \rQs - P\ r TsQ\ ~ P\TsU\ r 
+ P\Ts v *\ r Fs r + P\TsV* F 3rr (67) 


( ; u *xFlr + v 5x)Fsr + U* r F 3r + U*F 3n + V*,. 
= X F\ + v* lx F 3r + u* r F 3r 
+ u \F 3rr + v* r ^ 
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( 68 ) 



(“5«s* - «l«U + vIvL - vfvf* + JL r 5x - T lx )F 3r 
+ M*M */• - «l«T r + V 5 *V^ - vtvf r 
+ ^Z\Ts r -^T Xr = 0 (69) 

and 

2yl(T [x F 3r + T {r ){\+F] r ) + 2T x F ir F irr ] 

= l(7-D(Qi*G5 + Qi«to-Qiys*F 3r ) 

- (7-D2Qieix]/ 7 3r + [(7+l)(«?rQ5 

- v !r QsF-ir ~ Q5 v *F}rr + <2l«!r - Q\»5r^r 

~ Q\ vfcrr) ~ 2(7 - l)Ql(«lr - - vfo,,)] (70) 


The unknowns of the system are u\ x , v 4x , u* 4r , v* 4r , 
^5x> m 5 r> v 5x> v 5r> P4x> PSr> P%x> P4r> T 4x , T 4r , T$ x , 
T Sr , and the change in slopes of the slipstream and the 
resultant shock F irr . There are four shock equations ((67) 
to (70)) and four flow equations ((60) to (63)) in region 5, 
four expansion equations ((56) to (59)) and three flow 
equations ((60) to 63)) in region (4), and three equations 
((64) to (66)) valid at the slipstream h. The resulting 
system of linear equations is closed, with 18 unknowns 
and 18 equations. A matrix solution of the system yields 
the needed spatial derivatives to complete the solution for 
shock coalescence resulting in a shock, a contact surface, 
and an expansion wave. 

Concluding Remarks 

A method for finding spatial derivatives of flow 
quantities behind two-dimensional, nonuniform- 


expansion flow over a convex corner has been developed. 
Equations for the axial and radial derivatives of the two 
velocity components in the region behind the corner have 
been derived in terms of the curvature of the streamline 
behind the expansion fan. Taylor series expansions of 
flow quantities within the fan are used and boundary 
conditions are satisfied at the leading and terminal 
characteristics of the fan to obtain the first- and second- 
order solutions and the curvature of the characteristics of 
the fan. 

The method is applied to the case of sonic-boom ex- 
trapolation including the solution of shock coalescence 
with asymmetric effects. The calculation of the asym- 
metric effects requires finding the spatial derivatives 
behind the coalesced system. An outline is given for in- 
corporating the results derived in this paper into the 
system for those spatial derivatives. Application of this 
method will remove the restriction of continuous proper- 
ties through the expansion wave and extend the ap- 
plicability of the sonic-boom extrapolation technique to 
include expansion waves of finite strength. 

Langley Research Center 
National Aeronautics and Space Administration 
Hampton, VA 23665 
August 31, 1984 
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